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PACS 64 . 60 . F — Equilibrium properties near critical points, critical exponents 

Abstract. - We consider critical phenomena on heterogeneous small-world networks having a 
scale-free character but also arbitrary short-loops. After deriving the self-consistent equation 
for the order parameter and the critical surface, we prove that the critical behavior on complex 
networks is in fact infinitely robust with respect to the presence of arbitrary short-loops. 
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Introduction. — The simple but fundamental Curie- 
I Weiss equation is the prototype of all mean-field mod- 
[ els. For an Ising model with a coupling J /N defined 

• over the fully connected graph of N spins, the average 
\ magnetization m satisfies the celebrated equation m = 

• tanh{f3Jm + jSh), /3 = l/T and h being the inverse tem- 
[ perature and the external field. Similarly, Ising models de- 
1 fined on disordered graphs with a finite fixed (regular ran- 
[ dom graph) or random connectivity Poissonian distributed 
' (classical random graph), due to the tree-like nature of 
I the graphs, can be still exactly solved by using the Bethe- 

• Peierls (BP) approach [1], and the resulting average mag- 
I netization m obeys a more complicated mean field equa- 

■ tion which, however, presents the same critical behavior 
I of the Curie- Weiss equation with the classical critical in- 

■ dices; in particular m behaves as m ~ t^'^, r being the re- 
I duced temperature. In the last 10 years, a strong renewed 
1 interest towards mean-field models came from the study 

of cooperative phenomena defined on complex networks, 
i.e., random graphs having a heterogeneous distribution 
of links among the sites [2]. It turns out that, unlike the 
above random graphs, which are homogeneous, when the 
distribution of links is strongly heterogeneous, as happens 
when the degree distribution V{k) of the graph follows 
a power law Vik) ^ k~'^, with an exponent 7 < 5, the 
average magnetization m obeys an anomalous mean field 
equation giving rise to the critical behavior m ~ r^/h-^) 
for 7 > 3, or m ~ T-^/i^-i) for 3 > 7 > 2 [3]. As 
in the case of the homogeneous random graphs, also for 
the complex networks without short-range links, the tool 
to achieve the mean-field equation for m is basically the 



BP approach that is based on the tree-like assumption of 
the graph which, at least for 7 > 3, is locally true; in 
fact, if c is the average connectivity, the shortest loops 
are of length 0(log(A^)/ log(c)), so that in the thermody- 
namic limit there are no loops of finite length. The ana- 
lytical study of critical phenomena on complex networks 
has had a huge impact on the understanding of collective 
behaviors such as diffusion, cooperation and percolation, 
in biological, social or technological contexts. However, 
almost all these studies assume graphs locally tree-like (or 
at least hyper-graphs locally tree-like |3]) which in real- 
istic networks is almost never satisfied. Social networks, 
the neurons of the brain, the WWW and the Internet, are 
just a few examples in which the average clustering co- 
efficient C [2] is finite. More precisely, whereas networks 
having a hierarchical structure usually share a k degree- 
dependent clustering coefficient of the form C{k) ^ fc^", 
with a ~ 1, so that the most connected (and most im- 
portant) nodes are not clustered, there are other networks 
having C(fc) ^ 0(1) for almost any k. The former class 
includes e.g. some social networks, language networks, 
the WWW, and the Internet at the autonomous system 
level, whereas the latter class includes the Internet at the 
router level, the power grid, but also the brain. As dis- 
cussed in [5, the reason for this difference is related to 
the fact that in the second class wiring is expensive (eco- 
nomically or biologically) and the network, rather than 
hierarchically, is geographically organized. It is hence of 
fundamental importance to understand what is the role 
of a loopy topology in complex networks from the point 
of view of collective behavior. In particular: how does 
change the critical surface and - above all - are the an- 
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alytical results accumulated over ten years of research in 
complex networks robust with respect to the presence of 
loops? And if yes, to what extent? Notice that the crucial 
question concerns especially the loops of finite length but, 
due to the fact that the correlation length in these models 
remains finite, we are not allowed to neglect these short- 
loops neither near the critical point, so that the above 
questions are quite far from being trivial. 

At an intermediate level between loopy-like networks 
{e.g. finite-dimensional lattices) and tree-like networks, 
stay the so called small-world networks [B] in which we 
have an homogeneous distribution of links among the sites, 
but also a finite clustering coefficient, so that tree-like 
based techniques as the BP cannot be used. By using a 
completely different approach is however possible to solve 
exactly also these kind of homogeneous small- world mod- 
els at least in the paramagnetic phase (P) and to get their 
exact critical surface and behavior f7^. If Jq is the cou- 
pling associated to a given graph (£ojro), which in par- 
ticular may have short-range links and short loops, and J 
is the coupling associated to a number Nc/2 of uniformly 
spread long-range links (so that they form a classical ran- 
dom graph), then, for any c > 0, the mean- field equation 
in these models is given by 



def 



m = mo{p.jQ;ctm + /3h), t = tanh(/3J) 



(1) 



where mo^pjo; j3h) is defined as the average magnetization 
of the model in the absence of the long-range connections, 
with a short-range coupling Jq and in the presence of a 
generic external field h at temperature T = 1//3. Eq. 
([T]) is a natural generalization of the Curie- Weiss mean- 
field equation. It is easy to check that however, for any 
Jo > 0, the critical behavior of Eq. ([1]) is classical, regard- 
less of the local topology and clustering coefficient of the 
network (see also [8]). In this Letter we face the natural 
extension of Eq. ([T]) towards a large class of heterogeneous 
small-world networks that can be considered as zero or, at 
most, weakly degree-degree correlated. After deriving the 
equation for the order parameter and the critical surface 
(thermal or bond-percolative) , to which we devote also a 
couple of explicit examples, we prove that the critical be- 
havior (thermal or bond-percolative) on these networks is 
never affected by the presence of a local non tree-like struc- 
ture, provided the connectivity associated to such loopy 
structures be in turns non heterogeneous. 

Ising model on complex networks with loops. 

The family of models we shall consider is built as follows. 
Let {Lq,Tq) be any graph, £o and Fo being the set of ver- 
tices i = 1, . . . , A^ and links, respectively. Let us consider 
the Ising model defined on the graph {Cq,Tq) with a fixed 
coupling Jq and in the presence of an arbitrary external 
field {hi} 



Hf) — — Jn 



(ij)ero 



y^ hi<Ti 



(2) 



We will call this the pure model. Let us now consider the 
model obtained by removing randomly some links of the 
graph {Co, Fq) and by adding new links as follows. Let us 
indicate with Cg-ij = 0, 1 the adjacency matrix of the new 
graph in which some links of Fq have been removed. Given 
an ensemble C of random graphs c, c d C, whose links are 
determined by the adjacency matrix elements Cij = 0, 1, 
we define our heterogeneous small-world m,odel through 
the following Hamiltonian 



Hco,c — ^Jo 



(ij)6ro 



CO; 



IjUtUj 



y,- j^, 



•IjU lU J . 



(3) 



i<j 



The averages over the graph disorder (HI)-® will be in- 
dicated as ". The variables Ci_j specify whether a "long- 
range" link (and then a coupling J) between the sites i and 
j is present {ci,j = 1) or absent {ci,j = 0), whereas the vari- 
ables co;ij- specify whether a "short-range" link (i, j) G Fq 
has been removed (co;i,j — 0) or not (co;ij — 1). Both 
{co:ij} and {qj} are assumed to be independent random 
variables with distributions 



P0(C0;Jj) = (1 -P)'5c0;,.j,l +P5cO;i,j,0, 



(4) 



Pij{cij) = / {P'i,Vj) 5c,^,i + (1 - / {puPj))5c,^fl, (5) 

where p £ [0, 1], and the {pi > 0} are a set of normalized 
hidden variables, X)iP« = ^ 0110]: each proportional to 
the average degrees {ki} of the nodes £o in the absence 
of the links Fq. Usually, both /(•,•) and the {pi} depend 
on one, or more, continuous parameters ^ £l, and on A^. 
Given the mean degree of c, c > 0: 



E 



N' 



(6) 



we will assume that, for a continuous subset J7 C I, 
asymptotically in A^ 



I {Pi,Pj) = cNpiPj. 



(7) 



As a probability, Eq. ([7]) for f{pi,Pj) will be manifestly 
violated in X \ J^ whenever cNpiPj > 1. Note that, if 
Pi 7^ for any given N (a requirement which is true for any 
graph in which there are no isolated nodes), for A^ — > oo 
the terms cNpiPj tends either to or to cxi, therefore, in 
I\ J', the number of links {i,j) for which Eq. (O is not 
true for A^ large approaches 



A6 



def 



N = J2^ (cAp.pj- - 1) 



(8) 



where 9{x) = or 1 if x < or x > 0, respectively. If we 
parameterize A/jv as A/at ~ A^", and if the probability 'P{k) 
to have a vertex with degree k scales, for k large, as V{k) ^ 
k~^ , it is then easy to check that we have either a = 
for 7>3, orQ; = 3 — 7, for 7 < 3. As is intuitively clear 
(and as is possible to show rigorously [E]), the free energy 
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of ([3]), defined in the usual quenched way, has a leading 
OiN) term coming from the links satisfying Eq. ([7]), plus 
the 0{N°') rest. Therefore, in the thermodynamic limit, 
for 7 > 2, this rest becomes always negligible whereas, for 
N large but finite, the results we below present are exact 
up to corrections per spin which are near for 7 > 3, or 
0(iV2-7) for 2 < 7 < 3. 

As an example one can consider the choice 



The meaning of the order parameter m is quite natural: 



fiPi^Pj) = 1 



= 1 — e~^^PiPo 



Pi = 



(9) 



where jjl g [0,1). Eqs. © define the static model in- 
troduced in [9j. Note that the "fermionic" constraint 
that avoids to have multiple bonds in Eq. ^ is im- 
plicit. This constraint produces some weak dis-assortative 
degree-degree correlations for /i > 1/2 [12]. In the ther- 
modynamic limit N -^ cx), for fi G (0,1), Eqs. (J9]) lead 
to a number of long range connections per site distributed 
according to a power law with mean c and exponent 7 
given by 7 = 1 + l/^i, so that 7 e (2, 00). For fi e (0, 1/2) 
(7 > 3) the first of Eqs. ((H) takes the simpler form ^ 
while for /i e [1/2, 1) (2 < 7 < 3) it can be written 
as Eq. ^ only when i and j are sufficiently distant, 
ij ^ 7V2-i/ai^ while for lower distances, ij <^ JSf'^-^/t^^ 
we have Pij(cij = 1) ~ 1. 

The class of our heterogeneous small- world models given 
by Eqs. ©-([S]) with the only conditions c > and 7 > 2 
in I \ ^7 is very general. Note in particular that the graph 
(jCoiTo) is completely arbitrary and can contain closed 
paths of any length. We stress that the resulting net- 
work, union of the graph [Cq^T^) in which each link is 
removed with probability p, with the random graph c, can 
be seen as a gran-canonical heterogeneous generalization 
of the original small-world graph of Watts and Strogatz 
[5], though we here do not have a true rewiring. Since we 
let the probability p G [0, 1] and the mean c G (0, 00) arbi- 
trary, our way to build small- world networks is more gen- 
eral even for the homogeneous case pi = 1/A''. However, 
we can always restrict our class of small- world networks 
to the ones having a total average connectivity which does 
not change with p. To this aim we can choose c such that 
the number of links of the graph c is equal to the number 
of removed links of (£o,ro). Up to corrections 0(l/-\/iV) 
we can accomplish this for any sample by simply taking 
c = Cop, where co is the average connectivity of (>Co,ro)- 

Self consistent equation. — In the following, we will 
use the label to specify that we are referring to the pure 
model with Hamiltonian ^. Let moi(/3Jo, {/3^j}) be the 
stable magnetization of the spin i of the pure model ^ 
at inverse temperature /3. Then the order parameter m of 
the model ©-d?]) obeys 

m = ^mo.(;3jf ^ {A^PjCtm -f ph))p,, (10) 

i 

t'^^ tanh(;3J), 

io ''= tanh(/3J^^^) =-^ t&nh{pjQ){l~p). 



^Pi{<yi 



(11) 



((•) thermal average) while for any correlation function C, 
and for some positive 5, we have 

C = Co(/34^M^PjCtm-^/3/i})-t-o(7V"^), (12) 

where Co(/3Jo, {/3/ii}) is the correlation function of the 
pure model ^ at inverse temperature /3. The proof of 
(Unil-lini) is almost identical to the proof of Eq. dU [7], 
the only difference being related to the fact that, here, the 
probability pij (cy ) depends explicitly on the link (i, j). As 
Eq. dl]), also Eqs. ([T0l) - (fT2l) are exact in the P region and 
provide the exact critical surface and critical exponents, 
while m is correct up to 0{k/k'^) corrections out of the P 
region. A detailed proof will be reported elsewhere [11] . 

As an immediate consequence of Eq. PU)) we get the 
susceptibility x (= dm/d{/3h)) of the random model 



,(13) 



E. p. E, XO:i.j (/3 Jf ^ ; {Npictm + (3h} 
1 - ci^E.,j Xo:^J (l3J^^'^\{Npictm -f /3/i}) p,p, 



where xo;ij(/3>^o; {/3^j}) stands for the two-points con- 
nected correlation function of the pure model ^ at inverse 

temperature /3: Xo;i,i — {'^i'^j)o — {<^i)o{o'j)o- Note that, 
even if the model © with hi = h is translational invari- 
ant, the model ©-([S]) is no more translational invariant 
(except the case pi = l/N). Note in particular that x is 
not proportional to the sum of the connected correlation 
functions. In fact, from Eqs. (TTTI) and (rT2|) it follows 



XI P« [{<^^'^3) - (c^^X^j) 



(14) 



«J 



Critical surface (thermal and percolative). — 

From the self-consistent Eq. (flUl) we see that, in the 
thermodynamic limit, the critical surface /3c satisfies 



ciciV^XO; 



«J 



(/?,jf^0)p,p, = l. (15) 



For any given set ( Jq, J, p, {Pi}), Eq. (fTS]) gives the criti- 
cal surface of the model in the plane (/3, c). Furthermore, 
from Eq. (flUl) it is immediate to recognize that /3c < Pcq , 
with (3c — (3cQ only for Jq = 0, where /3^q is the critical 
temperature of the pure model ^ but with coupling Jq . 
The theory can be projected toward the limit /3 — >■ cx) 
where we get an effective percolation theory. Here the re- 
gion P corresponds to the region in which, in the thermo- 
dynamic limit, the parameters (c, co,p) are such that no 
giant connected component exists (m — 0). Note in par- 
ticular that, if Cqc is the percolation threshold of the pure 
graph {Co,Tq) (if cqc does not exist we can set formally 
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cqc = oo) in order to remain in such a region, the connec- 
tivity Cq = co(l— p) of the graph obtained from the graph 
{Co,^o) in which each hnk has been removed at random 
with probabihty p, must satisfy co(l — p) < cqc, otherwise 
a giant connected component aheady exists. From Eq. 
(fT5|) . for co(l — p) < coc, the equation for the percolation 
threshold Cc as a function of p is given by 

CcNj2xQ;r,j (tanh-i(l - p); 0) p,Pj = 1, (16) 



where we have used the fact that hm^_j.oo tanh(/3 Jq ') — 
tanh(l ~ p). Alternatively, Eq. ([T6| can be rewritten in 
terms of only graph elements as 

c,N J2 {k, + A/'ofi) P^PJ =1, co(l - p) < coc, (17) 

where Mq^^ a = 1 if, in the graph {Lq^Tq) from which each 
link has been removed at random with probability p, be- 
tween the vertex i and the vertex j there exists at least a 
path of links, and M^-lj = otherwise. 

Given p, if the condition co(l — p) < cqc is not sat- 
isfied, then a giant connected component is present and 
we can set Cc = 0. It is interesting to see in more de- 
tails the case in which we choose c = cqP so that, as we 
vary p, the total connectivity is fixed and equal to cq (as 
in a rewired small- world) . This study is important since 
it leads us to understand how the presence of short loops 
affects diffusion processes on general networks. In particu- 
lar, a strong interest regards the question: "in the presence 
of short loops how does the percolation threshold change?" 
If we set c = cop, from Eq. (fTT)) we get the percolation 
threshold cqc as a function of the rewiring parameter p 

4cVa^E ('5'^.■ +<i) P^P3 = 1' c'-iXi-p) < coc (18) 



From Eq. (|18p we see that p may have two effects on 
Cq\. . one tends to decrease Cq\. while the other tends to 
increase CqJ . In general CqJ decreases with p for p small, 
while grows for p larger. Let us consider the case in which 
Pi = 1; i.e., the classical small- world (no heterogeneity). 
In this case Eq. (IT5|) simplifies as 

c^^,Vxo(tanh-i(l -p);0) = 1, ci^^l - p) < c^,. (19) 

For example, if (£07 To) is the Erdos-Renyi graph (in the 
canonical representation), with mean connectivity co, by 
using in Eq. ^ xaifiJo^^) = [1 - co tanh(^Jo)]"^ (valid 
in the Po region) we get back obviously the well known 
percolation threshold c}^^ = 1, independently of the value 
of p. Depending on the problem, given co < coc, in general 
one can be more interested in reading Eq. (|17p either as an 
equation for p or for c. We will see later a simple example, 
while further analysis and applications to the percolation 
problem will be reported elsewhere [11] . 



Critical behavior. — We prove now that, if c > 
and (£0, To) is not a heterogeneous graph, then the critical 
behavior of the model ©-([71) does not depend either on 
(£o: To) or on (Jo, J,p, c). In [7], we have shown this result 
for the homogeneous small- world case pi = \/N and with 
p = 0. More precisely in such a case the critical behavior 
has been shown to be classical mean field for tQ>Q, while 
for ip < first-order phase transitions are also possible 
(see also [IS])- Here we will restrict the analysis only to 
the case tQ> Q. First of all from Eqs. ([T^ and ([15]) we 
observe immediately that the critical exponents for the 
susceptibility, above and below the critical temperature, 
are both equal to 1. Let us now turn to the analysis of the 
order parameter near the critical point. In the following 
we set ft, = and it is understood that we consider only 
the positive solution of (ITU|) tti > 0. For Jo = we have 

mo»(/3jf ^ {/3/ij}) - tanh(/3/i,) (20) 

so that Eq. ([IU|) strongly simplifies in 

m — g{m), g{m) = y ^ta.nh{NpiCtm)pi. (21) 

i 

The critical behavior of the model for Jq = 0, has been al- 
ready studied in [H]. Eq. ^ is equal to Eq. (21) of [Ti] . 
We recall that, if the distribution {pj} has a power-law 
character, we cannot derive the correct critical behavior 
of the system by simply expanding the sum g{m) term 
by term for small m. As shown in |14| . it is necessary 
to keep track of all the terms of the sum g{m). This is 
done by evaluating the sum with the corresponding inte- 
gral which gives rise to a singular term proportional to 
(m)^~'^ plus regular terms proportional to (to), (tti)'^ and 
so on, where 7 = 1 + l//i, if /i is the power law exponent 
associated to {pj}. As a consequence, when we solve Eq. 
(j^T]) to leading order in m, if T and r indicate the tempera- 
ture and the reduced temperature, respectively, we get the 
well known anomalous mean-field behavior: m ^ 0(t^'^) 
[i.e. classical mean-field) for 7 > 5, m ^ Oij^^'^'^^'^^) for 
3 < 7 < 5, and m - o(T-('^-2)/(3-7)) for 2 < 7 < 3. 
Note that the critical behavior of the order parameter 
m — ^j {<Ji)pi is different from the unweighted one de- 
fined as m = ^j (cri)/A^ when 2 < 7 < 3. In such a case, 
from Eq. (jl2[) one can use m ^ tm from which it follows 
that m - o(T-i/(3-7)) for 2 < 7 < 3. 

Let now be Jo > 0. It is clear that if the graph {Cq, Tq) 
is in turn a pure scale free graph with exponent 7', then 
the joined network will have an anomalous critical behav- 
ior characterized by the minimum between 7 and 7'. Less 
obvious is to understand what happens if (£o;ro) has 
a finite dimensional structure or some special topology 
with short loops. In particular we can pose the ques- 
tion: "does the critical behavior change by adding, via 
short loops, many paths between far spins, or may the 
critical exponent for m depend on Jq ?" Let us consider 
the self-consistent equation ([T0|) in general. The exact 
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expression of moi{j3jQ ,{l3hj}) for a generic non homo- 
geneous external field {hj} represents a formidable task. 
Note that, as above mentioned, to analyze the critical be- 
havior we cannot expand for small fields {/3hj}. However, 

since we have proved that, for Jq > 0, /3^o > f^c, in 
the region /3^q > /3 > /3c we can perform an high tem- 
perature expansion of m,Qi{j3JQ ; {/3/ij}) in the parameter 
ip = tanh(/3jQ ). It is then easy to see that 

+to [l-tanh2(/3/i,)] ^ tanh(/3/ij)4-0(ig), (22) 
jeJ\fo{i) 

where AfQ{i) is the set of the first neighbors of the vertex 
i in the graph (£o,ro). Plugging Eq. ^ into Eq. ^ 



m = g{m) + Ai(to) + A2(m), 



(23) 



where we have introduced 

Ai(m) =-^ioX!^* X! tanh(A^pjCtm), (24) 

and A2(to) is the rest coming from Eq. ^22\i which, near 
the critical point, |A2(to)| ^ Ai(to). Let us compare the 
bigger contribution Ai{m) with the "pure Jq = term" 
(7 (to) and let us focus on the simpler cases in which the 
graph {Co,^o) has a fixed connectivity |A/'o(i)| = Cq. Let 
us suppose first that cq = 1 (i.e., (£o:ro) is an ensemble 
of dimers). In general, given any normalized distribution 
Pi > 0, and any function f{x) > increasing with x, the 
following strict (for pi ^ 1/N) inequality holds 



E 



vtf{P]o{i)) <y2pdip^) 



(25) 



where jo(*) stands for the first single neighbor of i in 
{Co,Tq). Note that i — ^ joii) is a bijection on Co and 
that jo{i) ^ i. We can however formally enlarge the def- 
inition of jo{i) to include also the case jo{i) = i. The 
inequality (|25|) tells us that when we choose joii) = * we 
get an optimal overlap between the distribution {pi} and 
the function /(•). For the general case |A/'o(J)| = cq > 1, 
for any i, we can enumerate the cq neighbors of i in the 
same way as Jq (i), . . . ,Jq (i), so that for each upper 
index I (that represents an oriented direction) Jq ' {i) is a 
bijection on Cq on which we can apply (^5]) to get finally 






Co 

l — l i i 



By applying this equation to our case with f{x) — tanh(x) 
and for ig > we see that for any ?7i > 

< |A2(m)| < Ai(m) < tQCogim). (26) 

In conclusion from Eqs. (|10p and (E51) we get 

g{m,) < TO < (1 + i(,co) g{m) + 0{tl)g{m) + o{g{m)), 



where o(a;) is such that linij^^o o{x)/x — 0. In general, for 
to finite, it is possible to prove that |11) 

g{m) < TO < xo (/34^^ O) g{m) + o{g{m)), (27) 

where xo{PJq ^ 0) is the susceptibihty of the pure model 
^ with couphng Jq ^ and hi = 0. Since near the critical 

point, in the region /^^g > /3 > /3c, it is Xo(/3>^o ; 0) < co, 
we see that Eq. p7|) implies that the critical behavior of 
Eq. (jlOp remains always that corresponding to the term 
g{m), i.e. as if it were Jq = 0. Hence, g.e., if (£ojro) 
is a do-dimensional lattice and {pi} generates a scale-free 
graph with degree exponent 3 < 7 < 5, Eq. (P7)) tells 
us that near the critical temperature for the model ([3]) 
m ^ T^/h--^)^ regardless the values of Ja,J,p,c and do- 
We can also consider the case in which (£ojro) is a Pois- 
sonian graph with mean connectivity cq. To this aim we 
can start from the fully connected graph and remove from 
it randomly each of its N{N — l)/2 links with probabil- 
ity p = 1 — co/N; the resulting graph will be our Pois- 
sonian graph with mean connectivity cq. Since we have 
already proved that when (Cq^Tq) is the fully connected 
graph with a couplings 0(1/A^) the critical behavior re- 
mains equal to that of the model with Jq — (the effective 
coupling in this case being /3Jq = tanh(/3Jo)co/iV), we 
conclude that, also when {Co,To) is a Poissonian graph, 
the critical behavior on the total graph remains the same 
as if it were Jo = 0. 

Examples. — Concerning the critical surface, partic- 
ularly simple are the cases in which (Cq^Tq) is a set of 
disconnected finite clusters (points, dimers, triangles, . . . , 
or mixtures of them). Note that for such {Co,^o) there is 
no percolation threshold (or formally coc = c»). 

Let {Co, Po) be a set of N/2 disconnected dimers. This 
case represents the simplest example with Jo 7^ in which 
"^o^(/3Jo , {l^hj}) can be exactly calculated. We have 



tanh(/3ft,i) + to tanh(/3/ijp(i)) 



mo^[pJo : iPn.i) - ^^^^ tanh(/3/iO tanh(/3/i,„(,)) 



(28) 



where tg = tanh(/3 Jq ') and Jo{i) stands for the first neigh- 
bor of i. By derivation we get the correlation function 
Xo;i,j which, in the P region, takes the form 



XO;id{PJo ;0) = Sij +to (5»,jo(i). 

Therefore for the critical surface we have 



ctcN 



E^» +*oEp^P. 



%(*) 



1. 



(29) 



(30) 



With respect to the critical surface of the model with 
Jo = we see in Eq. ([50]) the presence of a new term 
proportional to Iq. How much this term affects tc de- 
pends on how the dimers are placed, i.e., on how we 
choose the first neighbors {jo(*)}- Since by definition the 
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dimers are not connected, in general for jo(*) we can take 
jo{i) — i + k,inodN where fc is a constant integer in the 
range [1, A^]- The exact evaluation of tc for N large re- 
mains simple only if k does not depend on TV or fc = 0{N). 
Under the choice © for the former case we have 



ci,(l + t,o)(l - m)'(1 - iV2^"i)(l - 2fi)- 



1, 



(31) 



whereas for the latter the term proportional to to in the 
thermodynamic limit becomes negligible. Of course, we 
re-find that for fi > 1/2 {i.e. 7 < 3) ic -5> 0+. 

Let (£0, Fo) be now a set of 7V/3 disconnected triangles. 
For simplicity here we assume fj, — 0, i.e., the graph dis- 
order is just the classical one (Erdos-Renyi graph) and 
we calculate the percolation threshold for the problem 
with fixed total connectivity that we have under the choice 
c = cqp, with Co = 2. For the susceptibility of the pure 
model with generic coupling /3 Jq we get 



xo(/3Jo;0) 



3g3/iJo 



-/9Jo 



jS^Jo ^ Se-z^-^o ■ 



(32) 



From Eq. (IT9t and (|32|) . by using the replacement /3Jo -^ 
tanh^ (1 — p), we get the equation for the percolation 
threshold Pc (which is solved for pc — 0.183406) 



o^3tanh (l—p) i p— tanh (1— p) 



(33) 



In the above examples we had do — 0. For a chain of N 
spins we can evaluate Tc from Eq. ([T5|) by using Xo;i,j — 
Vq . In Fig. (1) we plot simulations for the susceptibility 
X as a function of T for several system sizes and compare 
the location of the maximums with the theoretical T^.. 



Fig. 1: Plots of X as a function of T for a chain of N spins 
in networks generated via Eqs. ^ with c = 1, J = Jo = 1, 
/i = 1/3 (corresponding to 7 = 4), and p — 0. The vertical line 
comes from the solution of Eq. US]) with iV = 131072. 

Conclusions. — The results presented in this Letter 
are easily generalized to the case in which Jo and J are ran- 
dom couplings and in particular we find that also for the 
spin-glass critical behavior a similar robustness theorem 
applies. These, as well as other considerations concerning 



the existence of strong finite size effects (see Eq. (HU) will 
be reported elsewhere [TT]. Our analysis sheds also some 
light to the issues related to complex networks embedded 
in a metric space. Although we here do not give the long- 
range couplings J a metric structure, the graph (£0,1^0) 
is completely general so that it can be also a metric graph 
equipped with an arbitrary distance ||i — j|| , which is the 
natural distance through which the correlation functions 
via Eq. (1121) are expressed. An interesting generalization 
of the Ginzburg criterion for the relevance of critical fluc- 
tuations has been recently derived for Ising models built 
on annealed networks embedded in a metric space [TS] . We 
recall however that the Ginsburg criterion alone, when ful- 
filled, does not provide the critical exponent of the order 
parameter m, especially if a heterogeneous distribution of 
links, and then a continuous dependence on the parame- 
ters of this distribution, is considered. At the heart of our 
robustness theorem ([?7|) lies the general inequality (P5|) 
whose Ihs and rhs corresponds to the short-range and to 
the long-range couplings of the model, respectively. We 
see then that the existence of long-range couplings in a 
model makes its order parameter already optimized with 
respect to the addition of any possible regular, or weakly 
random (i.e. Poissonian like), structure; only the addition 
of further long-range links sufficiently heterogeneous may 
change the critical behavior of the system. This work was 
supported by SOCIALNETS. We thank A. L. Ferreira for 
many useful discussion. 
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